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The existing 256 chips long complex scrambling code ia WCDMA for high data rate users is 
not suitable for delay sensitive applications. A new fast code acquisition algorithm for 
WCDMA uplink is proposed. The complex part of the user specific short scrambling code is 
concatenated with leng^ 13 binary Barker sequence at the and 129^ chip. A Barker Code 
Match Filter is used at the Base Station for joint detection of multipalh signals of different 
Mobile Stations using Constant False Alarm Rate detection. For user specific code acquisition, 
a parallel search by Breadth-First Tree Search method using Truncated Sequential Probability 
Ratio test is used. The results obtained by simulation indicate that the mean acquisition time is 
much less than the full length searches. 
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Chapter 1 
Inttoduction 


1.1 Motivatton 

With the convergence of Computers and Communication systems there is a desire for greater 
bandwidth and more advanced systems while being wireless. The voice and data rates available 
in the first and second generation mobile communication systems no longer meet the growing 
demands in the twenty first century. Iligli speed latemet and multimedia services are the order 
of the day. The third-generation (3G) mobile systems (called IMT-2000) are betng designed to 
support wideband services, at data rates as hig^ as 2 Mbps, with the same quality as a fixed 
network [1] [2] [3] [4]. Everyone desires that the wireless systems should be able to provide the 
same QoS guarantees as that of the wired systems. In order to realiae this need, a new 
wideband wireless access technology incorporating recent technologes should be made 
possible. 

Wideband direct sequence code division multiple access (WCDMA) is a predominant 
technology for wireless systems of the third generation. WCDMA is desigied to offer ffeadible 
wideband services such as Internet access and real time QoS. However, these services are 
handicapped due to the physical limitations imposed by the wireless channel such as bandwidth 
constraints, multipath fading, noise, Doppler shift, multiple access interference etc. These pose 
a technical rballen^ in desi^iiag systems of higji bandwidth for mobile communication 
applications [5] [6] [7] [1] [8] [9] . 

To combat the effects of multipath fading, a Rake receiver is used to coherendy combine the 
mul%afb components in WCDMA systems. Theoretically, it may be possible to combine a 
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large nmnbet of paths imder no restdction of resources. This, however, is not feasible under 
most practical situations due to prohibitive complexity. In addition, certain real time 
applications are delay sensitive and demand daat the processing time taken be as low as 
possible. Hence, designing a low complexity receiver with minimum processing time is an open 
problem. 

1.2 Ptevious Work 

In the past few decades extensive research work has been carried out in the field of code 
acquisition and various methods have been su^ested by a number of them. Most of the 
researchers have considered AWGN channel modd with Multiple Access Interference (MAI) 
and Inter Path Interference ^I) modelled as white/coloured noise. However in advanced 
receivers the MAI and IPI are considered as stmetured and the correlation is taken into 
accounL limited work has been done in this fidd. 

In the case of parallel schemes, Rick and Milstein [10] present a parallel PN code acquidtion 
scheme for the reverse link of a cellular CX)MA system. Dong and Blostdn [11] sugg^ted a 
tree-search based searcher. The model of a terrestrial mobile commumcation channel 
incorporates the effect of MAI, shadowing, power control error, vehide speed, voice activity 
and sectorisation. The result shows that the system user capacity based on acquisition 
performance criteria is less than the capacity based on BSR criteria for certain ranges of system 
parameters. The effects of data modulation on the performance of parallel code acquisition 
scheme ate reported in [12]. 

In serial search method Deng [13] proposes a burst mode PN acquisition processor for non- 
coherent detection. This MF based PN code acquisition processor achieved a mean acquisition 
titnp of less than 10 bits over a range of chip — noise ratio (QNR) from 7 dB to 8 dB. Chang, 
Park and Lee [14] proposed a double dwell code acquisition system with a serial search that 



operates simultaneously ■with verification in long-code CDMA. For the search block, a MF is 
used for fast acquisition. For the verification block, an active correlator is used to reduce the 
complexity. The performance enhancement of this scheme mainly lies in reducing the false 
alarm probability by adding the verification stage while maintaining short acquisition time. 

Tree-Search detector has been studied as a low complexity approach in various research areas 
of wireless communication systems. In [15], a simple reduced tree search detection of the 
breadth- first type is applied to sub-optimal joint multiuser detection bit-synchronous CDMA 
systems over both Gaussian and tu^o-path Rayleigb fading channel by Wei and Rasmussen. Xie 
and Rushforth combined sub-optimum tree search algorithm with a recursive least-squares 
estimators of complex signal amplitude for joint Maximum Likelihood (ML) sequence 
detection and parameter estimation [16]. In [17], tree stmctured signal detection in V-BLAST 
system is proposed. It is found in [15], that tree search detector achieves near-optimum 
performance at a very low complesdty. Tree-search scheme in conjunction with sequential PN 
code acquisition to reduce prohibitive complexity of parallel PN code acquisition have been 
su^ested by Dong and Blostein [11]. However, the emphasis of their work has been more 
towards reduction of complexity and performance of gain combining technique when tree 
search is employed. Moreover, the system had inherent problems of rejection of large number 
of paths if null-hypothesis decision takes place at initial sta^s of the search. The paper has not 
made an attempt in studying the perfonnance of the detection technique and suggesting 
mediods for memory reduction and performance enhancement. There are short-cotrtings in 
the statistical model considered in their paper. 

1.3 Summary of Contributioiis 

The contdbutioas of dais thesis are: 

• Propose a new parallel code acquisition scheme with two different phases of operation. 

First multipath detection foHowed by user detection within each multipath- 
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® Propose change in the existing WCDMA uplink pilot channel structure by 
concatenating the imaginary part of the complex scrambling code with Barker 
Sequence for joint detection of multipath signals of all uplink users. 

® Use the Constant False Alarm Rate (CFAR) algonthm for joint detection of multipath 
signals of ail uplink users with the help of Barker Sequence. 

• Analyze and study the changes in autocorrelation and cross-correlation properties of 
the pilot signal as a consequence of the proposed change. 

• Propose a Tree search method using sequential detection for users present in each 
multipath. 

• Obtain analytical expression for the test statistics using Truncated Sequential 
Probability Ratio test being used in the Tree search. 

• Study the channel estimation being carried out for each user in a multipath during Tree 
search. 

• Analyze and study the mean acquisition time and avera^ sample size for the TSPRT. 

1.4 Thesis Outline 

In the foEowing chapters, we propose a new low-complexity receiver architecture for fast code 
acquisition in uplink WCDMA receiver in which path detection, user detection and channel 
estimation are also incorporated. The objective of this study is to achiev'e a low-complexity fast 
code acquisition in uplink WCDMA. 

We begb our study with a brief introduction to WCDMA communication systems in Chapter 
2. We will see the evisring WCDMA uplink system in light of the advantages and disadvantages 
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of code acqijisitioii. Following this we formulate the problem to be solved in this thesis work 
and the corresponding system environment. 

In Chapter 3, we have a treatment on the Constant False Alarm Rate (CFAR) method of 
detection. The chapter has an analysis on CFAR detection for Rayleigh faded signals and 
another analysis for PN sequence acquisition using CFAR detection for signals at the output of 
a correlator. Beginning with a Fixed Sample Size (FSS) test, we shall be following it by a 
treatment on Truncated Sequential Hypothesis Test (TSPRT) and its advantages over FSS. We 
win conclude with a description on Multistage Hypothesis Test (MHT), which is a modified 
form of TSPRT. 

Chapter 4 forms the core of the thesis work. We will suggest the changes required and see the 
resulting properties of the modified sequences. Next, the joint detection of multipaths by 
CFAR is covered. We then proceed further with an insight into partial correlation properties 
which win be instrumental for sequential detection. 

In chapter 5 we have the results of computer simulation and performance analysis of proposed 
algorithm. 

Chapter 6 summarizes the conclusion of this thesis and presents topics for future research. 
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Chapter 2 

Backgtoutid of WCDMA Uplink Receivet 

2.1 Introduction 

A spread spectrum CDMA scheme is one in which the transmitted signal is spread over a wide 
frequency band, much wider than the minimum bandwidth required to transmit the 
information which is being sent [9]. It employs a waveform that appears random to everyone 
other than the intended receiver of the transmitted waveform. For the convenience of 
generation and synchronization by the receiver, the waveform is pseudo-random, meaning that 
it can be generated by precise mathematical rules, but statistically it nearly satisfies the 
requirements of a truly random sequence. 

A CDMA system has several advantages over the conventional narrowband TDMA or FDMA 
systems; such as one-cell frequency reuse, inherent resistance to multipath fading, narrowband 
interference rejection, soft capacity, soft hand-off, simplification of channel alocation problem, 
ability to exploit silent periods in speecii voice activity, inherent message privacy and low 
probability of interception [18] [19] [20] [21] [22]. Though CDMA promises to be 
advantageous, there are many teciioical issues to be resolved. Synchronisation of PN sequences 
is difficult to achieve in the presence of multipath effects, Doppler shifts and fading. Fast 
power control is essential in the reverse link in order to overcome the near-far effects. Multiple 
Access Interference (MAI), Inter Path Interference ^I) and thermal noise are the capacity 
limiting factors in CDMA. The situation in uplink is further a^ravated due to the random 
offset between multipath signals of the same user and different users causing correlation 
between them. This emphasizes the requirement to desigp better sequences. This results in the 
MAI and IPI as not being Additive White Gaussian Noise (AWGN) but as a more structured 
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kiterfetetice thereby maimg the requirement of adTanced Multi-User Detection (MUD) 
techniques necessary. The advanced receivers provide better interference cancellation features 
resulting in better performance and network capacity improvement at the cost increased 
complexity of the receiver. Moreover, for packet data communication system, there is an 
additional requirement of fast code acquisition for real time Quality of Service (QoS) 
applications. Thus, there is a need to design a low complexity receiver airhitecture which is a 
trade-off between performance and complexity. 

2.2 Technical Features of WCDMA 

In WCDMj\ systems, there are sereral new reqturements: high spectral efficiency, high data 
rate, inter-system handover and multimedia co mmun ications with (Afferent QoS of multiplexed 
voice, data, image and video as well as different delay requirements from delay sensitive real- 
time traffic to delay flexible packet data. To meet these varied requirement, more sophisticated 
signal processing and communication techniques are proposed [4] [23] [24] [3]. These 
techniques indude: short complex scrambling codes concatenated by orthogonal variable- 
spreading factor (OVSF) channelisation codes [23] [!}, time multiplexed pilot channel structure 
and I/Q code-multiplexed pilot channel structure [25], fast transmit power control (WQ [1], 
pilot s 3 mibol aided (PSA) coherent dday lock loop (CDLL), PSA coherent code acquisition, 
coherent RAKE comb inin g, fast cell search under asynchronous operation, soft and softer 
handover, inter-frequency handover, blind bit rate detection, transmit divetdty, multiuser 
detection, turbo coding and smart antennas [26] [8]. 

2.2.1 Spreading Codes 

The spreading codes in CDMA may be dassified into short codes and long codes. A short 
spreading code has a period eqiKtl to the symbol duration and the long spreading codes have a 
period which is much longer than the symbol duration. In long code systems the corrdalion 
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between users changes from bit to bit and MAI is therefore random in time. This causes the 
performance for different users to be statistically identical Short codes have cross-correlations 
that remain unchanged over time and at times cause performance degradations over longer 
period of time. 



Figure 2,1: Spreading for uplink DPCCH, DPDCHs and HS-DPCCH 

A WCDMA system uses a combination of short and long codes and hence has multiple 
spreading. The data is first spread by an OVSF code followed by spreading by a PN sequence. 
This second level spreading is cel-specific in the downlink but common to al users of that cel 
and it is user-specific in the iq)lin]i. TTie short orthogonal codes (OVSF) ate caled as 
channelisation codes and the PN codes arc caled as scrambling code. Each transmission from 
either a BS or a MS is distinguished by die combination of a channelisation code and 
scr am b ling code. The PN codes in WCDMA are of either long or short type. The spreading 
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for uplink WCDMA is shown in Figure 2.1 [27] and the channelisation code is shown in Figure 

2 . 2 . 



Figure 2.2: Code-tree for generation of Orthogonal Variable Spreading Factor 
(OVSF) codes 


In Figure 2.2, the channelisation codes are uniquely described as where SF is the 

spreading factor of the code and k is the code number, 0 < yfe ^ SF-1 . Each level in the code 
tree defines channelisation codes of length SF, corresponding to a spreading factor of SF in 
Figure 2.2. 

The ^neration method for the channelisation code is defined as: 
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The leftmost value in each channelisation code word corresponds to the chip transmitted first 
in time. 

2.2.2 Sctambling Codes 

WCDMA uses two types of scrambling codes: the long codes and the short codes. The long 
scrambling sequences and Cj^^ are constructed from position wise modulo 2 sum of 

38400 chip segments of two binary ^-sequences generated by means of two generator 
polynomials of degree 25. Let x, and j be the two »-seqtjences respectively. The x sequence is 
constructed using the primitive (over GF(2)) polynomial + 1. The y sequence is 

constmcted using the polynomial X^^+X^+X^+X+1. The resulting sequences thus constitute 
segments of a set of Gold sequences. The sequence is a 16777232 chip shifted version of 
the sequence 

Since the primary concern of this thesis is relating to short scrambling codes for uplink 
WCDMA, we shall see in detail about this sequence. The short scrambling sequences 
and defined from a sequence of the family of penodically extended S(2) codes. 

Let « 23 ^ 22 - • -ttfl he the 24 bit binary representation of the code number n. 

The quaternary S(2) sequence 0 ^ ^ 16777215, is obtained by modulo 4 addition of 

tihree sequences, a quaternary sequence aQ) and two binary sequences and d(i), where the 
initial loading of the three sequences is determined from the code number n. The sequence ^(/) 
of length 255 is generated according to the following relation: 

- Zn(i) = aii) + 2h(/) + 2d(i) modulo 4, i = 0, 1 , . . ., 254; 

where the quaternary sequence 0 ( 2 ) is generated recursively by the polynomial 
+j:/4-2x‘+ 1 as: 

- fl(0) = 2«o + 1 modulo 4; 


a(i) = 2n/ modulo 4, i = 1 , 2, .. ., 7; 
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- aii) = 3a(i-3) + a(i-5) + 3fl(i-6) + 2a(i-7) + 3fl(i-8) modulo 4, f = 8, 9 254; 

and the binary sequence l>{i) is generated recursively by the polynomial 
as 

b(i) = n&n modulo 2, / = 0, 1 , . . 7, 

bii) = b{i-l) + b{i-J) + b(i-7) + fc(i-S) modulo 2, i = 8, 9, 254, 

and die binary sequence d(s) is ^crated recursively by the polynomial x^+x^+x^+x^+I 

as: 

d(0 = «i6+; modulo 2, j = 0, 1, 

dif) = d(i-l) + d{i-3) + d(i-4) + d(,i-8) modulo 2, £ = 8, 9, 254. 

The sequence is extended to length 256 chips by setting ^(255) = :^(0). 

The mapping from ^(/) to the teal-valued binary sequences and 0, 1, 

255 is defined in Table 2.1. 

Finally, the complex-valued short scrambling sequence is defined as: 

0) = mod256)t + j{- 1)' mod256)/ 2j)] 


where / = 0, 1, 2, . . . and U denotes rounding to nearest lower integer. 



Figure 23: Uplink short scrambling sequence generator for 255 chip sequence 
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An implementation of the short scrambling sequence generator for the 255 chip sequence to 
be extended by one chip is shown in Figure 2.3 


Table 2.1: Mapping from to Cshort,iA>i and '•= I* — > 255 


Zrih. 

Cstiort,1 

Cshott^A^ 

0 

-ft 

+1 

1 

-1 

+1 

2 

-1 

-1 

3 

+1 

-1 


2.2.3 Pilot Channel Structure 

The motivation for having a pilot channel is to have coherent RAKE combining, coherent 
code acquisition, channel tap weight estimation, power control and handoffs. There are two 
types of pilot channel in WCDMA system — common pilot channels and the dedicated pilot 
channels depending on whether the pilot symbols are shared by other users or not According 
to the structure of pilot and data symbols, it is classified as ttme-mulfiplexed pilot channel and 
parallel pilot channel structure. 

In parallel pilot channel structures, the pilot symbols ace transmitted with the data symbols at 
the same time. The parallel transmission can be T/Q multiplexed, code multiplexing or 
frequency multiplexing. Pilot channel shown in Figure 2.4 provides good tracking of channel 
status more so in fading channels. But the spectral efficiency is low and the pilot channel will 
cause interference to other users even when no data communications occur. In time- 
multiplexed pilot channel stmcture as shown in Rgure 2.5, the known Np pilot symbols arc 
periodically multiplexed into IVj, transmitted data symbols. Therefore, it is not suitable for 
delay sensitive real-time traffic such as voice and video transmission in circuit-switched 
systems. To allow for high speed data transmission over padret switched wireless networks, 
time- multiplexed pilot structure is ficvourable and it is proposed by Y Honda and M Mohar 
et al [28] [29]. The same shall be used in the thesis too. 
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Data Channel 
Pilot Channel 


Figure 2,4: Parallel pilot channel structure 


DPDC 


DPCC 



Figure 2.5 Time-multiplexed channel structure 


2.3 WCDMA Receiver Architecture 

2.3. 1 Basics of RAKE. Receiver 

In WCDMA systems, the transmitted signal band\»idth is typically lar^r than the coherent 
bandwidth of the channel. The hi^er chip rates in WCDMA communication systems give 
more multipath diversity than that in narrowband CDMA. Conventional modulation schemes 
require an equaliiier to undo the inter-symbol interference between adjacent symbols, CDMA 
spreading codes are designed to provide very low correlation between successive chips. If these 
multipath components are delayed in time by more than one chip duration, they appear like 
uncorrelated noise at a CDMA receiver and therefore equalhier is not needed. These 
uncorrelated paths are called resolvable paths because they could be separated by the 
correlation of a specific shift of the PN code [7]. 
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Figure 2.6: Uplink Rake Receiver for L path diversity 

RAKE receiver attempts to collect the time-shifted repEcas of the oti^al signal by providing 
separate correlation for each resolvable path and the output of each correlator is weighted to 
proxide a better estimate of the transmitted signal than is provided by a sm^e component. 
Demodulation and bit decision are then based on the weighted output of a bank of correlators. 
The RAKE receiver is essentially a diversity receiver designed for CDMA where the diversity is 
provided by the fact that the multipath components are practically uncorrelated from one 
another when their relative propagation delays exceed by more than one chip duration [7]. In a 
RAKE receiver there are correlators each of which are matched to different paths with 
different delays. In the event of a path facing fading, the other paths being uncorrdated 
undergo relatively less fading and the RAKE receiver can correct the decision values. This 
would not have been possible had there been only one corrdator. Decision is based on the 
combination of the separate decision statistics offered by the RAKE receiver. Thus, it provides 
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a forni of diversity which can overcome fading and thereby improve CDMA reception* The 
RAKE receiver uses different forms of combination for the decision statistics namely 
Maximum Ratio Combining (MRC), Equal Gain Combining (EGQ and Generalized Gain 
Combining (GGC). In order to realize a high speed implementation of a RAKE recei\^er, 
correlators are replaced by Code Matched Filter (CMF). With increase in data rates there will 
be more and more multipadis implying diat there will be a requirement of large number of 
CMF thereby increasing the complexity of the system. A typical RAKE receiver for an uplink 
in WCDMA is shown in Figure 2.6. The objective of the thesis is towards the design of a 
searcher for a RAKE receiver which continuously searches for changing multipaths and assigns 
fingers to the strongest multipaths. 

23.2 Channel Estimation 

Channel estimation problem in WCDMA receivers can be divided into two categories: 
multipath delay estimation and multipath tap amplitude and phase estimation. Theoretically, 
the channel estimation filter can be: 

• a predictor: if it uses only the past samples to estimate the current channel 

coefficients 

• a '‘filteF’: if it uses past and present samples for channel estimation 

• a smoother. if it uses past, current and future samples 

The removal of data modtilation can be accomplished in data-aided (DA), decision directed 
(DD) or non data-aided (NDA) manner [6]. In DA mode, the receiver knows some of the 
transmitted symbols which are called pilot symbols. Based on the received signal and known 
pilot symbols, the receiver can estimate the channel profile. The DD channel estimators utilize 
the decisions of die receiver to remove the effect of data modulation, where there is a decision 
feedback at the output of quanfizer. The NDA channel estimators (also called blind channel 
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estimatots) estiinates the channel ■without using data or decision. Instead, the statistical 
properties of the transmitted signals are exploited. 

Channel tap weights estimation can be based on linear interpolation, low-order Gaussian 
interpolation [30] [31] and Wiener filtering [32]. However, die Wiener filter requires the 
knowledge of statistical properties of the fading channel, i.e., fading covariance matrix, which 
is difficult to estimate in practice thou^ not impossible. Weighted Multi Slot Averaging 
(WMSA) ts generally used for channel estimation filter that simply averages the pilot symbols 
belonging to more than two consecutive data slots. 


2.3.2.1 Time-multiplexed Pilot Channel 

In die time multiplexed pilot chaimel, there ate pilot S 3 rmbols before Njj data symbols. 
Because fading remains almost constant over a period of Nj, pilot symbols, the simple 
averaging of consecutive pilot symbols improves the Sigpal to Noise and Interference Ratio 
(SNIR). The resultant channel estimate is given by: 




m=0 


( 2 . 1 ) 


which approximates the instantaneous channel gain at the time position 
n * + (Np — 1) * ^ in the t slot of the / path, where Tj (m, n) is the received symbol 

at data slot after RAKE combining. By averaging, the SNR of the channel estimate increases 
Np times fi:om the received SNR per symbol As in the case of slow fading, the channel gain 
remains static over a period of several slots, we can coherently add several channel estimates by 
a linear fiber with 2K taps to extend the observation interval Tlie channel gain varies slot-by- 
slot and hence the instantaneous channel estimates need to be wei^ted and summed. 
Therefore, the filter output is expressed as: 
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+ + ( 22 ) 

1=0 1=0 

where - is the real-valued weightmg factor and 2K is the observation inter\^al fepresented by 
die number of slots. Using the weighting factors ^s pre-stored in memory, the channel 
estimates (fi) are easily calculated from (2.1). In [25] the weighting factors optimised at the 

center of the slot are used for the reception of all the data symbols within the slot, that is jS ^ , 
is constant regardless of the symbol position of m. 

(m, n) = (n) m = , N^. + 1, N^, (2.3) 

By selecting the appropriate weighting factor accurate channel estimation is possible, 
particularly in slow fading environments, because a hxgc number of pilot symbols bdongtng to 
multiple slots can be used. If the fading is slow, the 2K —tap WSMA channel estimation filter 

^ . In the 

of fast fading, the channel gain changes even ■within the slot. Consequently, we change the 
wei^t factors symbol-by-symbol within the slot to better track fast fading. 

23.2.2 Parallel Pilot Channel 

For parallel pilot channel, pilot symbols and data symbols are transmitted in different channels 
at tlie same time. We also weight and sum the continuously received pilot symbols and the 
wd^t coefficients are updated symbol-by-symboL The filter output is given by: 

Wf (m, «) = MmT + nT^y, -1- hT) 

2Z 


increases the SNR of the channel estimate by a factor of 


/l: 


(2.4) 
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where is the real- valued wei^ting factor and 2Z is the observation interval represented by 
the number of symbols. 

The WMSA is a cascade of two filters: the first one is an integrator over a set of pilot symbols. 
It is followed by a linear filter with adaptive weights. Soft-decision statistics at the output of 
coherent RAKE combiner is represented at the symbol of it slot associated with the t 
propagation path. I = 0,1, , L — 1 as 


d (n, ^ (m, n) * iv, («) 

1=0 


(2.5) 


2.4 Problem Formulation 

Radio propagation in the land mobile channel is characteihied by multiple reflections, 
diffraction and attenuation of the sigial energy. These are caused by natural obstades/bodies 
such as hills, buildings, lamp posts etc resulting in multipath propagation. In the frequency 
selective multipath channel, the signal energy is smeared into a certain multipath delay profile, 
where multiple version of the transmitted signal are received with different amplitudes and 
times of arrival [7] [1]. 

In WCDMA systems, a wider bandwidth is allocated to support multimedia services and higher 
chip rates are used to multiplex the data streams of different users. An increased chip rate 
means reduced chip duration and path dday spreads over a largpr number of chip periods. 
Therefore, more resolvable multipath is available at the receiver, where multipatih components 
are spaced at least one chip period apart. At 3.84 Mcps, the waves travel 78 m in one chip 
duration (3*10®m/3.84Mcp5). Besides, multipath channel profile is dynamically time- 
varying which requires a RAKE receiver with a variable number of RAKE fing^ adapting to 
the time-varying channel state information (CSl). 
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An optimal EAKJE receiver, known as Full Search RAKE receiver, requires unlimited 
resources (CMFs or correlators) so as to combine aU of the possible multipath components 
[33]. Optimal performance is practically difficult to achieve because of prohibitive complexity 
and imperfect channel estimates. In the chapters to follow, we propose a new low complexity 
algorithm which causes a trade-off between optimality and complexity. The new algorithm is 
proposed for a WCDMA uplink receiver where we use a t\^^o-level detector. In the first stage 
the different multipaths of all users are detected and in the second stage detector the different 
users present m each multipath are detected and simultaneously channel estimation is done. 

2.5 System Model in Thesis 

In this diesis we consider an uplink WCDMA communication system working with Frequency 
Division Duplex (FDD). There are K active users in the system. Each user is presumed to have 
a maximum of L multipaths that may be combined in the RAKE receiver for multipath 
diversity. 


OVSF 

Codti Gan 



Code2 


Figure 2.7: Uplink WCDMA transmitter model 

The transmitter is simplified -with one Dedicated Physical Data Channel (DPDCH) and one 
Dedicated Physical Control Channd (DPCCH) with I/Q multiplexing. The simplified system 
is shown in Figure 2.7. Further since our thesis relates to code acquisition in the RAKE 
receiver searcher, we shall deal only with the base-band signal processing for DPCCH. 
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The following assumptions have been made: 

• Binary phase shift keying (BPSK) modulation is used at the transmitter with rectangular 
chip pulse shape. 

• Extended Gold code family with length of 256 PN chips is used as the short 
scrambling codes. 

• Coherent detection is used, aided by a separate pilot channel. 

• The multiple sigpal paths of each user are resolvable. An exponentially decaying 
multipath intensity profile (h'flP) is used. 

• Synchronous CDMA is achieved where the dominant paths of each i;user are 
synchroniiied in terms of symbol boundary. 

• Fractional portion of PN chip offsets of multipath propagation are ignored. 

• Fast transmit power control of the dominant path is used to compensate the path loss, 
shadowing and fast fading so that the near-far problem is avoided. 

• The number of resolvable multipaths is unknown at the receiver. 

• An isolated cell is considered with one base station in the center and K active users, 
therefore, no inter-cell interference exists. 

• Raylei^ fading channel with Additive White Gaussian Noise AWGN. The channel 


noise is modeled as AWGN. 



2.6 Uplink Signal and Channel Model 
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The model foi: the thesis and assumptiotis has been discussed in section 2,5. However, for 
simplicity, only the DPCCH with binary phase shift keying (BPSK) modidatioo and square 
pulse shape are assumed. The foUowing WCDMA system uses baseband discrete-time signal 
and system model. 

2.6.1 WCDMA Transmitter 

The transmitter uses direct sequence spread spectrum (DSSS) modulation. The information 
data symbols are spread by user specified short spreading PN code with symbol length and 
BPSK modulated. Finafly, it is transmitted on a radio carrier. The transmitted distal spread 
spectrum signal from user / in one symbol interval can be represented by: 

= ( 2 . 6 ) 

where the chip time index n — 1, G and G is the number of chips per symbol or 

processing gain, is the chip value of PN code of user i at time index n, € {—1,1} , 4 

is the transmitted symbol value at this symbol inter\’^al. The transmitted signal will be dispersive 
by the channel and it is subject to interference and AWGN. 

2.6.2 The Multipath Channel Model 

Frequency selective Rayleigh Fading is used to model the reverse channel from the mobile 
terminal to the base station. The midtipath radio channel for f mobile may be described as a 
wideband tapped delay line (WTDL) model with statistically independent time-variant tap 
weight : 


(2.7) 
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where L is the total number of midtipalhs, which is dynatnic because of the timewarytag 
characteristic of wideband wireless channels [6]. The index j represents one of JL multipaths 
experienced by the / user and is the relative integer chip time delay with respect to the first 
arriwg component for f path of / user. 

2.63 WCDMA Receiver 

In the multipoint- tO“point uplink CDMA channel, the received WCDMA signal is 
superposition of K channel outputs from K active users, corrupted by channel noise: 

K 

r(n) = ^ 5*^ (/i) \ iji) -h z{yi) (2,8) 

i=i 


where symbol * indicates convolution and ^n) is AWGN with a double sided power spectral 




density of - K is the total number of users. Hence, after RF analog downconvemion, chip 


pulse shaping filtering and chip-rate sampling, the equivalent received baseband discrete-time 
signal in one symbol period is obtained by inserting (2.7) into (2.8) yielding: 


Kn) = 

1=1 j=i 

where is the PN code of user i after a circular shift. Since code acquisition is aided by 

pilot symbols which are not data modulated, 4 ~ 

2.7 Code Acquisition of Multipath Signals 

In WCDMA system, code acquisition is \"ery important because reception of WCDMA signals 
is possible only after die receiver is able to synchronise die local PN code with the PN code in 
the transmission intended for it [34]. TTie other channel parameters may be estimated only 
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when the acquisition has taken place. Code acquisition is geneially earned out in two steps: 
acquisition and tracking. Code acquisition is the process by which the receiver attempts to align 
its local PN code s phase to the incoming signal within half or less than half of a-chip duration. 
Since the autocorrelation peak is much greater than the side lobes, a correlation between the 
local PN code and the received signal yields are estimate of whether the sequences are aligned. 

Many different types of PN code acquisition architectures exist, the most commonly used 
architectures are the active synchronous correlators and passive matched filter implementation 
[13]. For synchronous correlator implementation, each new offset value must be computed 
with integration over the symbol period which results in longpc acquisition time. For 
passive matched filter implementation, a code filter matched to the PN code sequence is used. 
The despread signal components appear sequentiafly at the matched filto: output, the passive 
matched filter can do the code correlations for different offsets much faster than the serial 
correlator even thougji output values are madiematically equivalent [35]. The disadvantage of 
matched filter is its high complexity. 

At the searcher in the RAKE receiver, the locally generated code sequence wil be phase 
matched with that of the received signal. Code acquisition is the process of deter m i n i n g the 
time delay of each path of each user. Strictly speakmg, determining the exact sampling time is 
an estimation problem. In order to reduce the problem to one of finite dimensionality, the 
delay is assumed to be a multiple of chip-time In other words, chip-synchronous 
environment is assumed. 

2.7.1 ConventioJial Matched Filter Approach 

In a conventional MF receiver, the received signal filtered by a pre-determined Code Matched 
Filter (CMF), where each possible delay of the PN sequence of the desired user is encoded as a 
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separate set of co-efficient. The received signal is match filtered by the PN sequence for the t 
user with integer time delay z^.of the / path and is denoted by Z-^ 

In section 2.2.2 we have ^ven a description on the generation of short spreading sequences or 
the complex scrambling codes of length 256 in WCDMA. These imaginar y part of the 
sequences are used for code acquisition and channel estimation. 

Let us consider the knagmary part of the complex sequence and represent it where 

f = — land n = 0,1,2,..., G — L Here Mis the number of sequences in the family (M 

_ 15777215) and G is the period of the sequence (G = 256). For a conventional IMF case, we 
have at the output of the CMF matched to user i with delay of j chips in a system with K users 
and each having L multipath where L is a random variable: 

Z, =|;r(n)c;?, 

n=l 


n=l k—l i=l 

kMJ 


( 2 . 10 ) 




where (2.9) is used and where 


G K 




M=I |fc=I f=| 

kM,i 




( 2 . 11 ) 


The signal from the f mobile -with shift will have a processing gain of G chips per bit, and 
the interference from other users can be modelled as AWGN type interference [6], The 
correlation statistics 2^ can be used to test the presence of die signal delayed by^T, seconds 
from the / mobile. In (2.11) the first term results in a combination of MAI (d i), IPI (if 
k = i and I ^ j) and the second term results in thermal noise. 
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Chapter 3 

Detection Techniques for Code Acquisition 

3.1 Introduction 

Signal detection pfoblems in wireless communications receiver design are generally M-aty 
hypothesis testing. There we wish to decide among M possible statistical situations describing 
the obsert^ations. For any given detection problems, there are a number of possible decision 
strategies or rules that can be applied. However we will choose a decision strategy that is 
optimum in some way for our algoridim. We shall employ two different strategies for the 
proposed algorithm. First, a Constant False Alarm Rate (CFAR) detection scheme is employed 
for detection of all possible existing multipaths. Second sequential detection scheme using 
TSPRT is used for user specific code acquisition within each detected multipath. 

3.2 Constant False Alarm Rate Detection 

There are two aspects of CFAR detection. First, we analyse the CFAR detection scheme for a 
Rayleigh faded signal. Second, die PN sequence acquisition using CFAR detection using a 
CMF. 


3.2.1 Conventional CFAR Detection 

The basic concept of a CFAR technique is diat the voltage of a test cell is compared to that of 
a set of reference cells. If the test cell voltage is “similar^’ to those of the reference cells, a 
signal-absent decision is made. If the test cdl voltage is not similar, a si^al-present decision is 
made. CFAR processor using range window for reference cells is shown in Figure 3.1. 
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Output 

Figure 3.1: CFAR Processor 

3.2.1.1 Methods for Obtaining a Set of Reference Cells 

There are several methods of obtaining a set of reference cells [36]. The proper method 
depends on the type of interference environment. One common technique^ shown in Figure 
3.1, uses a range interval bracketing the test cell as the reference cells. Use of this set implies 
that the characteristics of the interference do not change over the range interval of the 
reference and test cells. This assumption is true when the interference environment consists of 
the sum of thermal noise and constant-power jamming noise. The assumption is not true for 
environments that include ground clutter as the clutter cross section tends to change 
significandy in relatively small range increments. 



3,2.1.2 Detection Threshold of CFAR Processor 



Figure 3.2: Block diagram of the CFAR processor Figure 3.3: Signal detection 
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The received stgnal is square law detected and sampled in range by the AT-1 ran^ 
resolution cells as shown in Figure 3.2 The set of samples {x)i^ is processed resulting in a 
noise level estimate V. The estimate 12 is multiplied by a predetermined scale factor T residting 
in a pulse detection threshold Hjy 

3.2.1.3 Pulse Detection in CA CFAR Processor 

The sample from the test resolution cell Xo is compared with the threshold • The multipath 
signal detection is declared if the sample exceeds the threshold JTq (Fig. 3.2). 

H^=VT 

where T is the scale factor, maintained constant the false alarm probability, 
of die noise level in the cell xinder test 


( 3 . 1 ) 

V is the estimate 


N 

i=l 


where N is the number of reference cells. The output of the cell under test is compared to the 
adaptive threshold according to 


{Hf (multipath signal is present); if > /fj, 
(multipath signal absent); f 


(3.3) 


where H, and Hg are hypothesis that the presence or the absence of a target in the test celL 

3.2.L4 Fluctuating Signal 

Essentially most outdoor motale communication signals produce amplitudes that are Raylei^ 
distributed (exponential distribution of power; or cross section) [37] [38]. In general, the effect 
of fluctuation requires hi^«: sigpal-to-noise ratios for higri probability of detection and lower 
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\ alues for low probability of detection than those required witia no fluctuating signal Sweding 
[ 39 , 40] has considered four cases, which differ in the assumed rate of fluctuation and the 
assumed statistical distribution of the cross section. The two assumed rate are: 

® a relatively slow fluctuation of the signal, such that the values of for successive chips 

of a multipath are statistically independent but remain virtually constant from one pulse to 
tlie next, and 

• a relatively fast fluctuation, such that the values of are independent from pulse to 

pulse within one symbol duration (i.e., during the integration time). The first of the two 
assumed distribution for the received-signal voltage is of the Rayleigh form (3*4), which 
means that the mean amplitude of the signal has a probability density function ^ven 
by 


f 


] 




LSXpf- 


1 

— exp — 
V ^ y 






(3.4) 


where ^ is the average amplitude over all si^al fluctuadoii, s is the signal-to-noise ratio. The 


second assumed cross-sectioo density funchon is 




4jt 




— 2x 


4x 


^-2x^ 


-exri 


A) \ \ j 




H. 


(3.5) 


The first distributiofl (3.4) is observed when the received signal consists of many independent 
scattering elements of which no single one few predominate. Many mobile communication 
systems have approximately this characteristic at microwave frequencies, and seen largely in 
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complicated scattering environment, which ate usually of this nature |J>JLOS). The second 
distribution (3.5) corresponds to that of a path having one main scattering element that 
predominates together with many smaller independent scattering elements (LOS with 
scattering). In summary, the cases considered are as follows: 

• Case 1 Eq. (3.4), slow fluctuation. The chips received from the transmitter on any one 
symbol duration are of constant amplitude throug)iout the symbol duration, but 
uncorrelated from one symbol to another. 

• Case 2 Eq. (3.4), fast fluctuation. The pdf is as for Case 1, but fluctuations are taken to 
be independent from chip to chip. 

• Case 3 Eq. (3.5), slow fluctuation. The fluctuations are independent from scan to scan, 
but the pdf is given by (3.5). 

• Case 4 Eq. (3.5), fast fluctuation. The pdf is as for Case 3, but fluctuations are 
independent from chip to chip. 

• Case 5 Non fluctuating. 

The pdf for cases 1 and 2 is indicative of a signal with many scatters of equal amplitude. These 
are typical for complicated muldpaths like in dense urban areas. The pdf for cases 3 and 4 is 
indicative of a signal with one large scattered and many small scatters. 

3.2.1.5 Analysis of the CA CFAR Detectof 

According to rule (3.3) the probability of pulse detection may be found as [41] [42] 


0^ 


= Pi^>7V\H,)=]fyiv)dy]fix^lH,)dXa 

0 TV 
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where (f) is the pdf of the estimate V, and /(xjHj is the conditional pdf of t±ie test 
sample under hypothesis H,. 

The probability of false alarm is evaluated substituting 5 = 0 into expression (3.6). 

P,, = >TVlH,)='f/Jv)dvJ/Ujff„)^ 13.7) 

0 TV 

In other words, the probability of false alarm is determined in general by [43] 

P»=Ev{p(^«>TV\H,)} ( 3 . 8 ) 

wliich can also be written as 


FA ) 


exd 


V ^ 




\dx^ ^ = Eu 


exd 


^-7V^ 




■ = M,, 


^ 'J' \ 


(3.9) 


where denotes the moment generating function (MGF) of the random variable V. 

Similarly, tlie detection probability Pj^ is given by 

=£,{?(-»« >7V|H.)} (3.10) 

Since under the signal-present hypothesis id, the mean = .^ (l -I- 5 ) , we can determine Pp by 

simply replacing Aq with >^^(l-l-5) in (3.7). 


P„ =Mv 


( T '' 


(3.11) 


In conventional CA CFAR processor the noise level estimation V is formed as a sum of all 
outputs of the reference window. In this case the MGF of the estimate Fis defined to be [41]: 
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My{U) = M^{U) 


(3.12) 


where M ^{U) is the MGF of the random variable x, . The moment generation function may 
be found as Laplace transformation of the pdf f{x) 


M JjJ )= Jexp {-m)f{x)dx 


(3.13) 


The moment generation function of the random variable x, is: 






(3.14) 


The MGF of the estimate V (3.12) of CA CFAR detector is: 


MyiU) = 




(3.15) 


Assuming that the primary target in the test cell is fluctuating according to Swerling II model 
[39, 40], substituting (3.1 1) into equation (3.15) the probability of target signal detection is 


i + - 


(3.16) 


Setting i=0, in equation (3.16), the design expression for the probability of false alarm is 


PB=[i+r]- 


(3.17) 


The scaling factor Tis then computed from equation (3.17) thus. 


r=(PMr^-i 


(3.18) 
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As the number of reference range ceUs becomes large (Af oo), the probabiHty of detection 

approaches 


lim 




I + - 


1 + J 



(3.19) 


and the probability of false alarm approaches 

Pp,, = [l +r]“"' = exp(-T) ( 3 . 20 ) 

Equations (3.19) and (3.20) are the expressions which describe the performance of the ideal 
(tixed threshold) detector. Thus, for homogeneous clutter background environment, the CA 
ChAR detector is the optimum detector in the sense that its probability of detection 
approaches that of the ideal fixed threshold detector, as the number of reference noise samples 
becomes very large. 


The CA CFAR detector achieves ihe design probability of false alarm and a hi^ detection 
probability in a homogeneous background environment, that is, when the received noise 
samples are identically distributed and statistically independent. In a real environment however 
the noise samples may not be homogeneous. 

3.2.2 CFAR Acquisition of PN Sequence 

In this section we analyse the CFAR detection scheme applied to the pdf of the output of a 
CMF. An introductory analysis to CFAR has been given in [9] and in the above section where 
we have found the threshold for CFAR for a signal in multipath Raylei^ hiding channel 
Further detailed analysis has been given in [44]. The I-Q Non-coherent CMF and the 
Threshold Setting block for one particular phase can be ^ven as under. In tfa^ section we shall 
apply the threshold of CFAR in PN sequence acquisition. 



33 


The decision variable is square of a Gaussian Random Variable. It has a non-central 
distribution at the stage and central^' distribution at the stageH^ . 



Figure 3.4: CFAR Adaptive Threshold Detector 

To evaluate the probabilistic performance of a detector we need to model the output of the 
comparison filter at single sample time. Consider ihe received signal vector r to be a ALv? 
complex vector which corresponds to the last AT samples of the received signal Then 
accordingly 

r' =[r[n],r[n-l], ,r[n-(iV-l)]] (3.21) 

Let the filter vector, h, be the Nx/ complex vector containing the filter in corrdation form as 

h'=[M0],Ml], MN-m (3.22) 

The output of the filter at a single point in time is the inner product of the received signal 
vector and the filter vector y = h^x , where t rqjresents the operation of conjugate transpose. 

The filter vector is normalized to have unit energy |/tf = 1 , and can be one of N different filters 
of the system. 




(3.23) 
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We model the received signal as the combination of a signal component and AWGN, 


r = As + w 


(3.24) 


wliere s and w are (Nxl) complex column vectors. These are die signal vector ^ilot with data 
\alue of 1) and AWGN vectors respectively. The noise vector is complex multivariate GWN 

with ^^ero mean and correlation matrix — £{ ww ^ } = er^I. y4 is z real amplhudc dhat is 

Jot 

related to the SNR of the signal. The SNR of the received code is defined with respect to the 

peak response of a CMF, SNR = 1 2<J^ . For the peak response the signal vector, s, 
contains the code of the desired multipath (assuming that the pilot signal is unmodulated Le.; 
with a value of 1) 


s = b = M0],&[1],...,MA^-1]1 with = 1 

Since the code vector, b, has unit energy tihe signal vector s has variable eneigy between 0 and 
1 . The inner product of the signal vector, s, and the filter, h, define a location in a correlation 
function. The absolute value of the inner product can be defined as the correlation coefficient 

of the filter and the si^al vector, /) = jh^S^ = jpi? + jif^. The first part is the desired signal 

part and the next part is a combination of IPI and MAI. If i=J, then it is IPI, Le.; 
autocorrelation function and ifi ^ j , then it is cross-correlation between the pilot signal and 
other modulated data resdting in MAI. 

The constant false alarm rate is obtained by estimating the cneigy in the received signal and 
using this to adjust the threshold. Consider the case when no noise is present in the system. If 
we look at the ratio of P,,[n] over x[«] we obtain the CFAR normalized correlation fiinction. 


Comparing a location in this function to the fixed threshold parameterT, is equivaleat to 
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companng P^ln] to the adaptive &re^<MTrlnl The latki of wiU always be 

between 0 and 1 , regardless of die amplitude of Defining each value in these correlation 
functions as a CFAR correlation coefficient, Pj . Where Pj is defined as 

_ p' 

P } i~j' the offset from Ac peak response. (3.25) 

^ N 


Phe aim is that at high SNR/ no noise, the threshold should be so chosen that the peak 
response of the ACF will be detected while the range side lobes and cross-correlation values 
win not. As the energy in the noise increases a noise doud forms around each point in the 
correlation function. At low SNR, the uncertainty each location in the ACF becomes large. 
This will give locations whose Pj is less than a, a positive probability being detection, >0. 

Also this will decrease the probability of detecting the peak response, < 1 . We would wish 
to set a for a given and then evaluate P^ for varying SNR. The determination of and P^ 
consists of evaluating the probability of one random variable exceeding another. In the case of 
CFAR detector this is the probability thatP„[ni] > 7jc[n] . In the general case this is expressed 

as: 


P = 



>a 



(3.26) 


It is difficult to compute P because the LHS and RHS of the inequality are dependent (x). In 
order to overcome it, we will carryout a tonsformadon into a new co-ordinate system that will 
allow for separation of the problem into sums of independent random variables. The new co- 
ordinate system should include h as one of its member. The basis can be represented as matrix 
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# widi^Q — h* Since tiie colimins of ® fo rm a basis tibe matrix is unitary. 




The received signal can be represented as the linear combination of the basis set,r = ^a, 
where a is the (A/xl) column vector containing the coefficients of the basis vectors. The 
coefficients are obtained by projecting the received signal vector onto the basis set, a = <& V . 
The coefficients are independent Gaussian random variables with ^ new co-ordinate system) 
mean = and variance = <T^I . 


Substituting for r in the LI IS of the detection event procedures. 



(3.27) 


Doing the same for the RHS of the detection event gives. 


i- jr V| = — 

2 2 2 ‘ 


(3.28) 


When these representation are incorporated into the overall detection event the following 
simplification arises, 

> r^|rV||=i. I t = Tf(l-T) (3.29) 

The probability performance measures, P, now is the comparison of two independent random 
variables. 



(3.30) 
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Define lie nmdom vnmbie ^ wid,y = 0.,.., A^-l, Smce *e a, * a« independent 

Gaussian random variables the r- 's are independent non-central Chi-Squared random variables 
with tv^o degrees of freedom. The Moment Generating Functions (MGF) of the Vj 's is given 

by. 


hAu) = E\er'-^\- 


1 




exp 




-2 


1 + <T U 


(3.31) 


Fixpressing the LIIS and RHS of the detection event in terms of rj 's, we have 

N-\ 

X = v^andy = 5]v^. . The MGF of X is just the MGF of and the MGF of Fis just the 
1=1 

product of (he MGF's of the rj 's,J = 1,. . . JNF-l: 


hx («) = 


1 




exp^ 


ilAh"s|"« 


1+cr 


p.32) 


^y(«) = fl 


1 


Ut (1 4- (T^u) 


exp^ 


2 


A(^} 


2 

s| u 


\ + (T U 


(3.33) 


It is interesting to note that the MGF of X depends on osAy ,<T^ . The MGF of Y 

depends only on <y~ and the energy in the signal vector in the directions other than h. This 
means that the detection probabilities do not depend direcdy on the codes and filter of the 
system but only on the p , and thus the pj , that they produce. 
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By subtracting fl from both sides detection event can be expressed as comparing a Chi- 
squared random variable to a deterministic threshold, 

P = > tY] = P^X -tr > 0] p 34) 

This probability function is just one minus the cdf of a random variable. Since the MGF is 
known. Pis computed by solving the contour integration problem. 


c- ^ 


(3.35) 


Thus we obtain the probabdity of detection event as a function of a constant false alarm rate. 
We shall subsequently use these results for desigiing and analyzing our new algorithm. 

3.3 Truncated Sequential Probability Ratio Test 

In the previous section we have discussed the CFAR form of signal detection. In this section 
we will discuss another form of signal detection which will be used for user specific code 
acquisition within each multipadi. Among the three most common decision rules — Bayes, 
minimax and Ncyman-Pearson, we would choose Ncyman-Pearson [45] criteria because the 
least a-priori knowledge is required, in which a bound is placed on the false-alarm probability 
and the probability of detection is maximised within this constraint [46] [47]. TSPRT is a 
mixture of FSS and SPRT. Hence in the following sections, FSS hypothesis testing, sequential 
detection and MI IT will be analysed. 

3.3.1 Fixed Sample Size Test 

The basic physical observation model that we wish to conader is an observed waveform that 
consists of possible signal presence corrupted by additive noise and our objective is to decide 
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whether the sigaal is present or not by processing a collection of samples with observation 
window sii:e N taken from the observed waveform. 

Suppose the signal in (3.24) has constant amplitude and has a Gaussian distribution, we want to 
test a positive shift in the mean [48]. Let y = 3'oJ2’— realiisations of (i.d.d) random 

variables Y = FpYj, We denote 0 as the common mean of each of the ’s. Consider 

testing the hypothecs Hq, under which Y, has a pdf fiy-^o), against a shifted alternative H,; 
that is, consider the hypothesis pair. 

H,:Y,~f(y-a),e = e, 

versus 

Hr.Y,-ny-0\e = a,>a, 

for all 4 where / (y, )is the Gaussian probability density ftmclion of Y^, with variance . The 
Neyman-Pearson (NP) FSS test is obtained by testing N samples and the log likelihood 
function L^) [45]: 

> T=> chooseH^ 

< r chooseHf) 

where ^ is the observed realisation of the random variable 



Z. 


= ln 


/(y, -<9o) 


= {a,-aM-(a,+a,)f2)lv^ 


(3.36) 


Since we consider Gaussian distribution for the decision variable, hence we dete rmin e the first 
and second moments. The first moment of the Z- ^ven 0 is 
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f(h-0,) 


= jini 


./(}', -^o) 

= ie,-e,f{r-\)lv 


f (x - ff)dx 


(3.37) 


where r — {6 — 0 q)I{ 6^ — do). The second moment is 


mg=E{Zf \0) 




= J 


= In 


I 


/U-^o) 


/ (x - 0)dx 




V 


V 


and the condittotial variance is 




r^.-^oY 


V » 


= <T" 


(3.38) 


(3.39) 


To reduce subscdpts we define jUy = /ig^ and //<, = . Then we have 

fi, = -{e,-e,fliv-=-ii, (3.40) 

The sample size iV and the threshold r are pre-chosen [48] so that the test has error 
probabilities P(choosmg j Mq true) and P(choosing true) of a and 1-^ determine the 

threshold, r, 

T = {«r) (1 - m(W(Mo - A)) (3-41) 

and number of test samples. 




(3.42) 
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whete 'PC*) represents the standard (nortnalbzed) Gaussian distribution function and 
^-1 

(®) is its inverse [48]. 

33.2 Sequential Detection 

We will now see a different type of hypotliesis where we will cany out hypothesis testing over 
each chip and truncate die test once we achieve the required level of confidence. Unlike the 
FSS test, there is an alternative approach where the desired performance is fixed and the 
number of samples is allowed to vary in order to achieve this performance. Namely for some 
realizations of die observation sequence we may be able to make a decision after only a few 
samples, whereas for some other realizations we may wish to continue sampling to make a 
better decision. A detection scheme that uses a random number of samples depending on the 
observation sequence is generally known as sequential detection. Wald’s [49] sequential 
probability ratio test is obtained by testing, at the nth sample. 

< jB Hq P(Zj » Z2 

>A 

i=l 

e (J?, A) => take another sample, (3.43) 

where the boundaries ^4 and B are chosen so that the error probabflities are AT and 1-fi. The 
sample size M = min{» : y" , z, € iB,A)} is now a random variable, and average sample 
number depends on the actual distribution of Z; , / = 1, 2, 3,. . which depends on the value of 


0. This can be represented diagrammaricaJly as shown in Figure 3.5. 
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Figure 3.5: Depiction of Sequential Detection 

The maifi benefits of using the sequential probability ratio test ^PR*T) lie in die following three 
facts: 


• Sequential detection achieves same probability of false alarm and detection as a fixed 
length test 

• The average test length for the sequential detection is smaller than for a fixed length 
test 

• The reduced average test length can be exploited to reduce computation complexity. 

An obvious disadvantage of the SPRT is that occasional long tests may result if the observed 
data are ambiguous. As the Wald-Wolfowite theorem [46] implies the average test length of 
SPRT is much shorter than FSS. However, these occasional long runs may not be practical for 
many applications. A practical compromise is a truncated SPRT (T-SPRT) where a finite-stage 
truncation exists [48]. As a trade-off, the optimality property in terms of minimiiang average 
sample sb:e (ASS) is lost by such a truncation. Nevertheless, the increase in ASS is usually 

sligh t. In the truncation stage, n = M , test 
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>t* 


{5.44} 


3.3.3 Multistage Hypothesis Test 


In the section on Sequential Detection we have seen the T-SPRT was expressed in terms of the 
log likelihood function, which is impractical for realistic system design. In terms of the received 
samples, tibe T-SPRT can be rewritten as [47]: 


r 


ty. 




< choose and terniinate 


> \ choose and terookiate 
else contmue to stage n+l 


(3.45) 


in stage «, 1 < n ^ iV, . At the iV, stage, we test 

n 

ly. 


i=l 


2 : Af, => chooseH^ 
else chooseHf^ 


(3.46) 


where the and JB„are thresholds at the stage n and N/is the truncation stage. In the 
following, we ^ve the transformation step by step. 


n n 


2=1 2=1 

vs 

A 

= E<t-(«.+e.)/2) 

i=l 

vs 

Av- 

(«, -»,) 

III 

vs 

Av^ 

(«,-«.) 

III 

vs 

Av' 

(«, - So) 


+ 1 - 


( 3 . 47 ) 


Here we define the upper threshold at stage n as: 
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Av- +<9o 


Similarly, the lower threshold at stage « can be expressed as: 
" 2 


where 1 < n < and the tnmcation threshold at stage can be expressed as: 




(3.48) 


(3.49) 


(3.50) 


From equations (3.48), (3.49) and (3.50), die 2N, —1 parameters for an IV^-stage test are 
reduced to four parameters: A, B, truncation point r and truncation stage The T-SPRT can 
be viewed as a mixture of SPRT and FSS tests. If Cg and r, are mixing constants each on [0,1], 
using Wald’s inequalities [46, 48], we set 


A = ln 

l-(l-c,)P« 

. (I-C„)P;. 

II 

r (i-c,)Pm ■ 

Ll-(l-Co)Pp. 


(3.51) 


(3.52) 


iV, =[»?-' (coP,) + 'F-'(c,P„)f 


M^-Mo 


and 


(3.53) 





^.54) 


where Pf, = CX, and P^ = I — . 


The design may be optimised by varying I’g and ly. values of £g and q near 0 yields a test amitar 
to SPRT. Alternatively, values of Cg and q near 1 may result in similarity to an FSS test. In the 
test design, <t, jUg, and ju, are known; Pp and ace nominally chosen design values of the error 

probabilities. Once Cg and q are chosen, test is completely specified by A, B, Nj and r. The 
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design guarantees tihiat for independent and identically distributed samples; die actual 
probability of false alarm and missed detection will be lower than nominal values. The 
performance of the MHT will be canicd out subsequendy. We shall analyse the probalhlity of 
MHT reaches a certain stage and the probability of accepting either hypothesis. 

3.3.4 MHT fof Rayleigh Fading Channel 

We shah now extend MHT for reception of Rayld^ faded signals at the output of a correlator. 
By Central Limit theorem, die pdf of the output of the correlator, when there are number of 
users wih have a Gaussian distribution. Moreover, if the signal has undeigone Rayleigh fading, 
then the joint distribution at the output of the correlator be a combination of Raylei^ and 
Gaussian distribution. 

The pdf of the amplitude of the received waveform is Rayleigh distributed. The chip energy E^- 
is multiplied by where g is Raylei^ disttibuted. Hence, 

= (3.55) 

where <7^ = £[5'^ ] and ^['j is the expectation operator. Let T = f ^ then we have 

1 f-r) 

\^R J 

where /^(jOis the pdf for Rayidgh disttibution, (Tg is the variance of the distribution. The 
average chip energy becomes 

(3-57) 

At the output of the correlator, we now have the following situations: 

• Under Hg, hypothesis, there are two possibilities: 

o there is no signal, the pdf is aero-mean Gaussian 

o there is Rayldi^ faded signal that is uncorrdated at the output of the cocedator 
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• Utidef H, hypothesis, the Rayleigh faded signal is present and it is correlated at die 
output of die correlator. 

Under null hypothesis and in the first case, it is a pure noise-only path that causes die 
thresholds to he higher "while in the second case, the Rayld^ faded signal is present but 
uncoirelated at the output of the correlator causing the thresholds to be lower. As will be seen 
later, in our system we shall use two level of detecdon and the MHT will be used in the second 
stage detector pro'tided the first detector (CFAR) decides that a signal is present. Moreover, in 
the uplink, there will be some or the other baseband signal present and hence we reject the no 
signal case. For the analysis of MHT in the case of Rayleigh fading, we consider only the 
second case under null hypothesis. 

Then as ^ven in [5^ the pdf considered in section 3.3.1 for the Gaussian case will now he 
modified by wd^iing the si^al power or equivaleody 6^ by y. Then we have 


Ss 


Hi -Y, ---^expl 




r.y\ 


f { y -}^\0 = Oi>0o 


The NP test variable will now be ^ven by 


Z, = In 


fiyi-0o) 


= ln 
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(if) 
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—exp 
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exi^ 

U ff 
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2V^ JJ 


-JLf/f -A /y 
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(3.59) 
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The difference between the Gaussian NP decision variable in (3.36) and the joint Rayleigh and 
Gaussian variable in (3.59) is the multiplication by the factor y. Accordtn^y, (3.47) wiH now be 
modified as: 


n 


J .z. 


1=1 



■vs 


Av^ + 


(?.6G) 


The upper threshold for stage » gets modified as: 


(3.61) 


and the lower direshold gets modified as: 


S - 

■ M-e.) 2 

The truncation stage threshold is gjven by: 


(3.62) 


(3.63) 


Next, we analyse the changes occurring to the mean, mean square and variance of the decision 
variable Z; gjven in (3.37), (3.38) and (3.39) respectively. The first momeit of Z- for a ^ven dis: 

\'\ 




-I'la niffv 




+^ o )' 


( 3 . 64 ) 






where r = {0 — J^o)/ ~ ^o) • second moment is ^"en by: 



48 


m,^E\^]j9]=E 


/ 

=4rW i'"-i;m+9o)+ 

U I 




r(^i +^'o) 




= ■^(3, -«„)’(«’ +«")+^(«, -9„) 




/ 






^{6, -e^f^^id, -e,f +^o)+ 


W V 

f ^2 , r 


7^ (^} + ^^) 
4 


2A 


u 


tW-So)’ +^W 


=4w -«'.)= +4w-'’.)'f'-7 

t> u V 2 


a2 


(3.65) 


The conditional variance is modified firom (3.39) and ^ven by; 




-^.)= -^.)*('-f) -(^w -^.){'-|) 

_ f 


Y 

y 






(3.66) 


The modified values of JH^ and (7^ from(3.64) and (3.66) can be used in (3.53) and (3.54) for 
obtaining the values of iV, and T . 
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Cliaptet 4 

Proposed Algorithm for Code Acquisition 

4.1 Introduction 

Parallel code acquisition schemes are computahonaHy intensive. Here we propose a new 
algorithm for parallel code acquisition in uplink WCDMA which has lower computational 
complexity. A searcher is employed by the Rake receiver to detect different multipaths in order 
to provide multipath diversity for the traffic channels. The BS provides the timing adjustment 
information to the MS [1]. This ensures that the multipath signals from different MS arrive at 
the BS at dose intervals. We first carryout non-coherent joint detection of multipath signals at 
the BS using CFAR technique. For this we concatenate the existing scrambling code of 
different MS with length 13 Barker sequence. TTie Barker sequence concatenated signals of 
different MS with same relative delay are jointly detected at the BS by a CMF matched to the 
Barker sequence. Then in each detected delay-path, we carryout user-specific coherent code 
acquisition by tree search. 

4.2 Tree Search Algorithm 

The tiumbct of CMFs for parallel search employing brute force approach would be number of 
users multiplied by die number of possible phases for each user. Instead we propose an 
alternate method of code acquisitton. This algorithm uses the short spreading codes of existittg 
users to construct a binary tree and then the tree is pruned breadth fkst for code acquisition 
using TSPRT. This has lower computahonal complexity since TSPRT requires much lesser 
number of samples than a FSS test [49]. 
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4.2.1 Code Tree Constructioii 

Instead of testing all possible circular shifts of all users PN codes independently by a bank of 
correlators/ CMFs, we propose to test them jointly and utilize common testing statistics. Some 
shifts of PN codes have initial chips in common; simultaneous hypothesis testing is possible by 
organizing all users PN codes into a code tree. This is achieved by first enumerating all possible 
PN codes widi all possible circular phase shifts of existing users in a code book. Then we 
proceed fhrou^ the codebook, codeword by codeword. Each codeword is processed chip by 
chip: each chip value is encoded in a code tree node. 

In a WCDMA uplink system, let us consider K users existing. Let codebook C denote a matrix 
with K rows. Each row corresponds to a user in die system with lengdi iV chips where N is the 
length of test For the convenience of algorithm description, we identify each node at a 
particular depth level of the tree search by the rule: 

• Node (/j denotes the/* tree node at tree depth level L 

• Node (Jy ffs left child is labeled as node (fkl, Z/l). 

• Node (iy jfs right child is labeled as node (/Fl, 2^. 

• / = 1 JSJ, \<j< 2' , where N is the number of chips over which the TSPRT test is 

carried out 

The sequential detection test can be done independently. 

The procedure of constructing the code tree works recuisivdy. The chip values of codebook C 
are labeled row by row according to the labeling rule. For row i, chip /. 


If CJr] = -K then generate a left child 
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• If CJ/] [4 = 1, Aen generate a right child 

The last row of the codebook consists of leaf nodes of the tree, i.c., no children are 
g^erated at the end of each row. The depth of the tree is, therefore, IV level 

The problem with this system of code acquisidon is that at the early stages of the sequential 
test if a rejection occuis, then a large part of the tree is not searched and this leads to higher 
missed detection. At the same time if a selection occurs at an eady stage, it will lead to a higher 
false alarm. We will surest an alternative method of tree search which will reduce both the 
miss detection and false alarms and also limit the size of the codebook. In order to limit the 
false alarm rate and miss detection, we wil restrict the size of the codebook. This is achieved 
by constructing the codebook only for existing users with zero circular phase shifts. We then 
carry out a two level search — first we have non-coherent detection of different multipaths that 
are existent in the uplink system and then we carryout coherent tree search within each delay 
path. For this we concatenate the short spreading codes of all useis with Barker sequences 
which have got excellent autocorrelation properties and can resolve multipath. As per the 
existing WCDMA standards [1], in the uplink the BS gjves timing adjustment command to the 
MS to adjust its transmit time. The timing resolution should be within 1.5 T, [IJ. This implies 
that multipaths of different users with comparable relative delays will arrive at the BS at close 
intervals and therefore make code acquisition by joint detection feasible. 

4.3 Joint Multipath Detection 

Barker sequences have excellent autocorrelation properties that can be utilized for multipath 
detection. The same is explained in the next section. 


sp© A 
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4 . 3.1 Barker Sequences and its Properties 

We define Barker Sequence as a sequence of +rs and -I’s of some finite length k such 
that the correlation function C(t) defined by [50] 

CW = satisfies |C(r)| < 1 for 7?^ 0 . (4.1) 

la (4.3) we have ^vea the expression for the aperiodic autocorrelation function. We now 
define the periodic or the even autocorrelation function (EAC) as: 

^(r) = C(r) + C*(iV,-r) (4.2) 

and the odd autocorrelatioa function (OAC) as: 

m = C{t)-C\N,--r) (4.3) 

where is the length of the Barker Sequence. The OAC function affects the output of the 
Correlator when the information symbols change over one integration interval, while the EAC 
function affects die output when the information symbols do not change. Thus, when the 
binary information symbols are equiprobable, both the EAC and OAC functions are equally 
impotlant in the system design and performance analysis. In addition to the above two 
properties there are two other interesting properties. Namely, the Merit Factor — which 
specifies the ratio of energy of the ACF main lobe to the energy in the ACF side lobes 



and die HNV which specifies the ratio of the ACF value at m = 0 to the ACF out-of-phase 
value with msudmum absolute value 
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HNV = V/K 5* OmodiVg (4,5) 

= max {C’(^)} • (4.6) 

There are a total of six binary Barker sequeoces. For our algo rithm we shall use leng th 13 
Barker sequence which will cater for a delay spread of 13 chip-time. The sequences and its 
correlation properties is as shown in the Table 4.1 and Figure 4.1. This is very important from 
our point of view since Barker sequences have a very good merit factor and their HNV equals 
the length N. This g^ves a measure that by how much will the in-phase components is above 
the out-of-phase components. The hipest known merit factor of 14.08 is that of length 13 
sequence. The correlation properties of binary Barker sequence of lengtli 13 are as shown in 
Figure 4.1. 



Figure 4.1: Correlation properties of Bariser Sequence of length 13. 
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Bgifeer Se^aueacc; atiid 'Propcrd.e$ 

MF 


Seq: {+1, +1, +1, +1, +1, -1, -1, +1, +1, -1, +1, -1, +1} 


13 

ACF; (13, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1) 

14.08 


PCF: (13, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1) 



OCF:(13, -1, 1, -1, 1, -1, 1, 4, 1, 4, 1, 4, 1) 



Table 4.1 : Correlation properties of binary Barker Sequences 


4-3.2 Concatenation of Short Spreading Codes by Barker Sequence 


In the previous section we have seen the excellent correlation properties possessed by the 
Barker sequences. These sequences can be used to detect multipalhs over smaller number of 
chips. When the receirred signal is passed through a Barker correlator the output at the zero 
time-shift is 13 dB higher than other shifts of a length 13 Barker sequence. Tlius, it will be 
possible to detect the different multipath signals with the help of Barker sequences. 

Let lis denote the length 13 Barker sequence by the vector 


ilXlXl-l-UX-iX-U) 


(4.7) 


for ti= U 13. The {cf } 
modified and represented as: 

fcf }= (fn, ) 


in section 2.7.1, which is the scrambling sequence, will be 
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1 < < 13 

where 14<n^ <128 . 

142<tij <256 

The transmitted digital spread spectrum signal from user / in one symbol interval can now be 
represented by: 

s,(n) = c,!% («) 

Ihe received baseband discrete-time signal in one symbol period is obtained from (2.10), (4.8) 


and (4.9) as: 
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K L 


rin)^ 

/=! M 


(4,10) 


This received baseband signal is passed through a CMF matched to Barker sequence with delay 
of p chips. The discretc-time output of the Barker CMF in one symbol period when M users 
have paths with delay of p chips is ^ven by: 


® A 

2, “E'-Wi..,. 

11=1 

»=1 is=l Jjs} 

^ B ^ 

M=l 1=1 

Jjj^ ^ 

+ +\-128-r, ) Az + 
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kl^ip 
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•ilk’ + C(T)iJ(n-T,)+C(r)<J(n-128-ir,))/ff« +z(n)^„_,^ ) 


i=l 


Jk = li=l 

kMp 


= E (4 f - T ^)- h 3 in-m - T ^ » 

i=l 

+ E +C(T)Sin-Ti) + CiT)Sin-i2^-Ti)]fi,^ + zin)b„_,^ ) 


fe=l 1=1 


M 

= £ ( Ap ) ■»■ - j 28 -T, )) + iA?, 


-Rfhcre IN p =E E((^r + ) (^-12) 


k=l f=l 


is the resulting interference and 
rip r«/-% 


is the cross-correlation between Barker sequence and the scombling code. The first term in 
(4.1 1) is due to the exact match resulting from a total of M users with the same time delay of p 
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chips* The second term is due to the balance of K-M users with path delays not matched to the 
CMF. 

Comments: The following are observed: 

• Tlie first term will be 13 dB higher than die second term with the C(t) terms and has 
been seen in the properties of Barker sequence. 

• The cross-correlation term in (4.13) can have a normalised value in the range fO, 1]. 

This means that even in the case when there is no exact match, diis teixa could lead to 
situations of false detection. Thus, there is a need to reduce this cross-correlation term. 

4.3.3 Decoffelated and Concatenated Sequence 

In tibis section we will address the problem of reduction of cross-correlation as given in (4.13). 
We shall now carryout modification on the concatenated sequence given in (4.8). 

In order to reduce the cross-correlation betw^een Barker sequence and the scrambling 
sequence, we will decorrelate the scrambling sequence part andc^|'). This is achieved by 

decorrelating the oti^al scrambling sequence at the transmitter with a filter that satisfies the 
following criteria: 



{d„ 7 = (0.076923, 0.076923, 0.076923, 0.076923, 0.076923, - 0.076923, 

- 0.076923, 0.076923, 0.076923, - 0.076923, 0.076923, (4.14) 

-0.076923,0.076923) 



The de-correlated and concatenated sequence can now be represented as: 


}= fc, > > K-m ’ * a„ ) 


( 415 ) 


where * represents convolution. The equations (4.9) and (4.10) will now be given by: 
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r(n) = + z(n) 




For the same given conditions as in section 4.3.2, we now have: 

n=l 

«=j i=i j=] 

Vs M 

= 1,1. K Wn -T,)+S(n-m- J-,)) 
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(4.17) 
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(4.18) 



where IN ^ = S E (ip¥^ + CiT)Sin -T,) + CiT)S(n - 128 - )) + zin)b„_,^ ) ( 4 . 1 9 ) 
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Figure 4.2: Decorrelated concatenated scrambling sequence 


In (420), we have an interesting obser\^ation. The cross-correlation between the scrambling 
sequence and Barker sequence is reduced mainly to the correlation between Barker sequence 
and the decorrelator filter. This is shown in Figure 4.2. The maximum value of the correlation 
is unity when the relative delay between Barker sequence and decorrelator filter is zero and 
otherwise the value is either zero or 0,0769. In either of the conditions the value of p^^^will 

depend on the summation of the original scrambling sequence over length 13. Due to the 
balance property of PN sequences, this summation will be a small value. When the signal at the 
output of the Barker CMF is passed throu^ a square law detector, the value of this cross- 
correlation will be very small and merge with the background noise. Thus, this modified 
decorrelated and concatenated sequence will have a very low probability of false alarm rate. 


In section 3.2 we have seen two separate cases of CFAR detection: first case is of a signal in 
Raylei^ fading channel and the second was for PN code acquisition. In the case of uplink 






WCDMA, the CFAR detection scheme for midtipath detection will be a combination of the 
two cases. After having reduced the correlation of out-of-phase signals to mere background 
noise, we can assume, without loss of generality that the detection scheme w^iH be only of 
former kind and neglect the latter kind. The CFAR detection scheme followed will be that of 
Figure 3.4 where the governing equation is obtained from (3.18), (4.17) and (4.18) 

z z 

j|z/ >(('’» r*’ - (4.21) 


Cross-correlation of length 13 Binary Barker sequence and deorrelator filter 



Figure 4.3: Cross-correlatioii of length 13 binary Barker sequence and decorrelator 
filter 

Comment: The concatenated and de-correlated sequences have been obtained by 

modifying the existing scrambling sequences so as to enable detection of multipath signals in 
uplink WCDMA with low probability of false alarm. Since the ori^al scrambling sequences 
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have been transformed/ modified, dieir autocorrelation and cross-correlation properties should 
be studied and compared with that of the original sequences. 

4A User Specific Code Acquisition 

In (4.15) we have obtained a transformed sequence that will be used for code acquisition in our 
algorithm. After a multipath has been detected by (4.21), we commence user specific code 
acquisition within each multipath by Tree search using MI IT. For this we use the scrambling 
sequence from the 14^ to 128*^ chip. The code acquisition procedure will be carried out twice 
in one symbol duration of 256 chips. This will be done as: 

• carryout multipath detection bet^^^een 1^ and 13^’ chip and 129^ and 141*^ chips. 

• once a multipath is detected at the end of 13^' chip, then carryout detection of users 
present in each multipath from the 14^' to 128^^ chip. 

• if a multipath is detected at the end of chip, then carryout detection of users 
present in each multipath from the 142^ to 256^ chip. 

The scrambling sequences prior to concatenation were binary and subsequent to concatenation 
and decorrelation are no longer binary. The consequence is that the codebook will become 
very complex and will be multi-valued. In the next section we wiB see how the sequences can 
be restored back to binary value at the BS. 

4.4.1 Inverse Filtering 

The WCDMA baseband signal processing is linear. The CMFs at the transmitter and the 
receiver can be designed using tapped delay line (TDL) filters. The weights of the Barker TDL 
filter and decorrdator are obtained by inverting the vectors g^ven in (4.7) and (4.14). We shall 
use the same analogy in designing an Inverse Filter for reverting the decorrdated part of the 
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scrambling sequence to binary values. This will enable us to carryout tree search using a binary 
tree. The inverse filter to be designed should be inverse of the sequence in (4.14). Let us 
represent the Discrete Fourier transform of the vector in (4.14) by D(k). Let the inverse of the 
transfer function be represented by (P(k)y\ The two transfer function will be related as: 
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Figure 4.4: Proposed code acquisition algorithm 


D{k){Dik)y=i 
{D{k))-^ = 


D(k) 


(4.22) 


The time domain representation of the inverse filter is given by; 


/ l_ -,)_r0.86667, 0.86667, -1.3,0.86667, -1.3, 0.86667, 0.86667, 

¥„i-F {[Dik)} j-l -1.3,-1.3,0.86667,0.86667,0.86667,0.86667 


(4.23) 


At the BS the chips arrive serially. The incoming chips are first passed through the CFAJl 
detector. If a hit occurs indicating the existence of a multipath, the CFAR outputs a logical 
signal that enables a block that wiU carryout die tree search. The received baseband will pass 
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through the Inverse filter before being input to the tree searcher block as in Figure 4.4. This 
will restore the part of scrambling sequence back to the original signal corrupted by MAI, IPI 

and noise. 

4.4.2 Partial Correlation of Scrambling Sequence 

Tree search involves MHT wherein a hypothesis testing is done after each sample and the test 
is terminated when the test statistic crosses either the upper (selection) or lower (rejection) 
threshold. As against this, there is the conventional match filter approach of section 2.7.1 that 
does FSS test. For FSS type of tests, the test statistics will depend on the correlation of PN 
sequence over its entire length while MHT will depend on partial correlation. In this section we 
shall analyse the partial correlation of the original scrambling sequences. 

In section 3.3.3 and 3.3.4 analysis of MHT is gven. In this section we shall derive bounds for 
partial correlation only for the tree search part of the sequence and not the entire modified 
sequence. The partial correlation of a PN sequence is ^ven by [51]: 

BJkJ) = ■£ 0 < *,/ < -1 (4.24) 

i=l 

which is a function of the initial phases and 4 the length of the correlation w and and jb are 
the two different user codes. m— p, then it is autocorrelation else it g^ves cross-correlation. 
When m = p and k- I, then is true and R^{k,l) is linear. When k^l, under Hg, then we 

have R^(k,l) in the worst case situation, which is approximated as linear. We obtaia by 
simulation the approximate expressions for bounds on the mean and variance of R^(k,l). 

W 

1=1 



Variance of Partial Correlation > Mean of Partial Correlation- 


(4.25) 
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Figure 4,5: Mean of partial correlation of scrambling code 
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Figure 4.6: Variance uf partikii cuneiatioo of scrambling code 
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, ^ W W 

mi c%c% .IlKTamoaisgivenby 

1=1 js=l 

0.0328>v^ 

V =h^0.012741 — J (4.26) 

These values will be used for MHT variables and for desigoing die thresholds for 
tree detection. 

4 A3 Data Stnictiire and Algorithm Description 

In this Stage of code acquisition we will use MHT as described earlier to prune the code tree. In 
case no users are detected in a multipath implies that it is a noise only path and a fake alarm 
has been raised. Because multiple hypothesis are tested jointly^ some users that are deemed to 
contain noise-only components will be elimnated simultaneously while others that are deemed 
to contain signal energy are detected simultaneously. The lower T-SPRT thresholds B, and 
upper thresholds and truncation stage threshold are designed using equations (3.61), (3.62) 
and (3.63) by setting the parameters as in (4.25) and (4J26). 

Because we expect to complete code acquisition in one symbol period, the maximum number 
of stages in the MHT is less than the number of chips per symbol JV, therefore we rewrite: 

N, =min(N,,N) (4.27) 

where is calculated by equation (3.53), (3.65) and (3.66). The following information is stored 
within each tree node: 

• The current test stage / Le., die tree depth level of node j. 

• PN code value of tree node J, 

• The test statistic of its father node where i is the tree depth level and J is the tree node 


identifier. 
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• A linked list that contains all paths that intersect this node. 

The tree search part of code acquisition is a combination of tree search, code correlation and 
sequential detection. We denote, Z,. ^as the T-SPRT test statistic of tree node (i, j), r(/) as the 

signal chip sample to be processed and j is the input test statistic of the father node of tree 
node (qj). The algorithm is executed as follows; 

1 . Initialise test statistic value at the tree root node (tree depth level 0) as Zg j = 0 and input 
this value to tree nodes (1 ,1), (1:,2), Le. ZFj j = ZFj 2 = 0 . 

2. Recursively execute the test in depth level / by breadth-first tree search. All the tree nodes 
at tree depth level / are tested- The test statistic value at tree node (L, j) is computed by: 

2u=C,,*r(i)-iZF,, (4.28) 

• IfZj j < Sj , the subtree from this tree node « discarded from further processing. 

• If Z. j > A; , recorded the detected paths intersecting with this tree node. 

• l(Af<Zfj<Bi, output Z, j to its child nodes at stage / + 1 and proceed to the next 
tree depth level, i.e., set / = / + 1 and repeat. 

3. If the final stage is reached, Le., i = N,, and if all the tree nodes in this stage are tested, 
tlien terminate the test. 

4.4.4 Chamiel State Information 

If the output of the CMF/corrdator used for a fuU length search yields a hit then the 
expectation of the output of the comelator given in (2.10) under the condition Zy > T ^es: 

e{z„}=-L(A(G + W,;) 



The second term is ne^gible as compared to the first term in the event of a bit occurs. 
Therefore, expectation the output of the CMF/correlator jields the channel estimate. In the 
case of a sequential search, the correlatioo can be considered to be “an incremental chip-sized 
code match filtenn^’. This means that the output of the CMF is checked at each time index 
and if it crosses the upper or lower threshold, a detection or rejection occurs. If it does not 
cross then the output is checked at the next lime index. Hence, an expectation of the 
expression in (4.28) results in the same value as in (4.9) when the decision statistics is greater 
than the upper threshold Given that ^ > .4,. in (4.28) we have: 


rk,}=-(Cj*'-0)+zF,,) 

I 


(4.30) 


which is the channel estimate 


for the ^ path of the user. Thus the output of fhe tree search 


averaged over the stage at which it detection gives the channel estimate. 
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Chaptet 5 


Peffomiatice of Code Acquisidoii Algorithm 

5.1 Introduction 

In tile previous chapters we have developed the necessary theoretical background for the 
modified code acquisition scheme for WCDMA uplink. It is difficult to evaluate the 
performance of the sequential tests because their calculation rapidly becomes unmanageable 
even under certain simple conditions [52] [51]. For example, when the sequential test is used in 
a noncoherent receiver with a Gaussian noise channel, it is difficult to analyse due to the 
unavoidable modified Bessel function in the likelihood function [53] [51]. The problem is 
further compounded by the fact that very litde literature is available for the study of uplink 
code acquisition under structured interference. U(&ik code acquisition is rather complicated 
even under AWGN channel [9]. Thus, the only other alternative available for studying the 
perfomiance is with the help of computer simulations. 

5.2 Software Experiments 

In order to get a realistic result for the uplink WCDMA code acquisition, we have designed the 
experiments as per the WCDMA specification TS 25 series [54]. We have taken a systems 
approach rather than a pure mathematical approach. The simulations have been earned out on 
Madab 6.5 and Simulink 5.0. 

5.2.1 Uplink Transmitter 

The u plink transmitter design is as shown in Figure 5.1. The transmitter consists of a pilot 
sequence generator which consists of a trial of continuous ones at 15 Kbps. This is then spread 
by the scrambling code at the rate of 3.84 Meps. The scrambling code type is different for 
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different expetimeiits (as has been described in the previous chapter) and varies &om user to 
user. The signal is next scaled by a factor § as has been shown in Figure 2.1. This baseband 
signal is then passed through a transmitter pulse shaping filter with a root-raised cosine (RRQ 
roll-of- factor of 0.22, oversample factor of 8, filter sampling frequency of l/(3.84e-6) Hz and 
has 96 filter taps. After this the signal passes throu^ lie multipath c hann el. The multipath 
channel uses a fading profile with a maximum Doppler shift of 40 Hz (pedestrian case) [54] to 
fade the signals. Each user has two or three muMpaths. As has been mentioned eariier, without 
loss of generality, we assume that timing adjustment and power control are existing in the 
system. Therefore, the users have the first path either with 0 or 1 chip delay, second multipath 
of 2 or 3 chip delay and third multipath of 4 chip delay. Thereafter, the signal passes through 
an AWGN channel which adds noise to the si^aL 



Figure 5.1: Uplink transmitter used in simulation 
5.2.2 Uplink Receiver 

The WCDMA uplink baseband receiver used for the simulatioo is as shown in Figure 5.2. The 
multipath faded and noise corrupted signal from different users is added at the recdver firont 
end. The received signal is then passed through a receive pulse shaping filter with the same 
parameters as the transmitter pulse shapang filter. The upsampled signal is then downsampled 
to 3.84 Mcps. Thereafter, the receiver signal processing is as shown in Figure 4.4. 
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Figure 5.2: Uplink receiver baseband signal processing used in simulation 

5.3 Peffofmance of Modified Sequence 

In section 4.3.2 we have analysed die concatenation of short scrambling codes that are 
described in section 2.2.2 widi Barker sequence of length 13. In this section we will evaluate 
the consequence of this concatenation by simulation. 

5.3.1 Concatenation by Barker Sequence 

In order to resolve multipath signals at the BS, we have proposed concatenation of the 
scrambling sequence with Bariser sequence. The sequence is as ^ven in (4.10). A wilful 
correlation is introduced between the signals of different users to enable joint detection. For 
this part of the experiment we have a sin^e user system with three multipaths under Raylei^ 
fading channel with path delays of , IT^. and 4T^. , where T, is the chip time and has relative 
mean power of -3dB, -6dB and -9dB. The output at the Barker CMF is ^ven in Figure 5.3. 

Comment: We observe that due to the cross-correlation between the 14” — 128* and 142"^ — 
256* chip and the Baiker CMF results in false peaks. This is as gjven by (4.13). This will cause 
false alarms in the system. These false peaks need to be removed. 
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Output of Barker Correlator at Reciver for pure sequence 



Figure 5.3: Cross-correlation of spreading sequence and Barker sequence 

5.3.2 Performance of Decorrelated and Concatenated Sequence 

We now investigate the sequence that has been obtained by further modification and given by 
(4.17). The scrambling code is modified by decorrelating the 14* — 128* and 142"“^ — 256''’ chip 
of the scrambling code. The simulation is carried out with same parameters as in section 5.3. 1 . 
The squared output at the Barker CMF is shown in Figure 5.4. The statistics of the simple 
concatenated sequence given by (4.8) and concatenated and decorrelated sequence of (4.15) at 
the output of Barker CMF is given in Table 5.1. 

Comment: In the previous section we have noticed that the concatenated sequence gives 

a large number of false peaks. In Figure 5.4, we observe that as given in (4.20) the false peaks 
have been largely reduced resulting in lower false alarm rate. A comparison of the statistical 
properties of the two sequences substantiates the fact that decorrelated and concatenated 
sequence will perform better in CFAR detectors. 
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properties of the two sequences substantiates the fact that decorrelated and concatenated 
sequence will perform better in CFAR detectors. 



Figure 5.4; Cross-correlation of de-convolved spreadir^ sequence 
smd Bsuker sequence 


Statistical Properties 

Pure Sequence 

Decorrelated Sequence 

Mean value 

12.1493 

0.9307 

Median 

9 


Standard deviation 

1 16.1744 ^ 

1.4636 i 


Table 5.1: Comparison of statistical properties of modified sequences. 


5.3.3 Correlation Properties 

In the previous section we have seen the benefits of decorrelated and concatenated sequences 
when it comes to CFAR detection at the output of a CMF matched to Bariter sequence. Ihe 
decorrelation so carried out over the short scrambling codes should be checked for any 
imdesirable effects due to tiiis transformation. The short scrambling code provides ideathy to 
each transmitter in the up link and their correlation properties should be retained. Any 
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imdesirable change in the coixelatioii properties can cause interference to the other channels. 


Autsocorreiatlon properties of original and modified sequence 

AAittteorr«f<ttian otf eomiktex part {Mira saq AutocorrfatfoR of compMx part of concotenatetf & docwrralatect seq 





Figure 5.5: Autocorrelation properties of pure short spreading sequence and 
decorre^ted concatenated sequence 

We shall now see the complex, real and imaginary autocorrelation and cross-correlation 
properties of pure sequence (without deconeklion and concatenation) and the modified 
sequence. 
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Cross-cfsnielation properties of original and modHled seciuences 
Cross-correlation or con^kfex part of pure seq Cross-correlation of complex part of concatenated & decorreiated aeq 



Cross-correlation of real part of pure seq Cross-coirelatlon of real part of corHudensled & rtocoiielaled seq 



Shift-— > 


Cross-corrolaUon of Imaginary part of puro seq Cross-coneladon of ctmtplax part of concatenated & decorreiated seq 



Figure 5.6: Cross-correlation properties of pure short spreading sequence smd de- 
convolved concatenated sequence 

In Figure 5.5 we observe that the autocorrelation properties of the pure and modified 
sequences are maintained except for the ima^ary case (which has been modifiec^. The two 
other short peaks are due to the Barker sequence and are not of much sigoificance during tree 
search as this part of the sequence is not used in tree search. In Figure 5.6, we observe that 
there is an increase in the cross-coccdatioa in the c^e of the ima^nary parts due to the 
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presence of the same Barker code concatenation for aJl the sequences in the family. Under 
similar reasoning as above, we can ignore the marginal excess correlation. 


5.3.4 Spectial Properties 



Table 5.2: Comparison of stabstical properties of power spe^rum 


the spectrum is noise like in both the csaes. The statistical properties of the power spectrum 
show very little changes. We, therefore, conclude that the properties of the modified sequence 
do not vary much and the added benefit is that we get better performance of multipath 
detection using Barker sequence. 
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5.4 Multipath Detection Performance using CFAR 

In this section 'we will analyse the perfotmance of the CFAR detection system. First we analyse 
the system for a sin^e user case with varying SNR and false alarm (threshold for CFAR) and 
then we study its performance for multiuser case under Raylei^ fading. The former test is a 
rather deterministic test to ascertain performance of the receiver design. The latter test is for 
the actual WCDMA code acquisition. 

5.4.1 CFAR Detector in AWGN Channel: Single User 

We have seen the properties of the decorrelated and concatenated sequence in the preidous 
section. The probability of miss detection for various false alarm rates is simulated for different 


Plot of CFAR Probabili^ of miss detection for Barker Sequence in AWGN 



Figure 5.8: Probability of miss detection versus SNR for different false alarm 

values of SNR in an AWGN channel The threshold for the different false alarm rate is gjven 
by (3.18). Using this value as the threshold for the CFAR detection, we first study the 
performance of the detector in an AWGN channel for a single user case. The probability of 
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detectioii is g^ven by (3.16). The results of sunuladon for probability of detection versus SNR 
for different values of false alarm are gven m Figure 5.8. The plots in the figures indicate that 
with increase in SNR, the probability of miss detection decreases and the probability of 
detection increases. For lower values of false alarm rate, there is lower probability of detection 
as more path detections are missed owing to higher threshold value for detection. 



Figure 5.9; Probability of detection versus SNR for different false alarm 
5.4.2 CFAR Detection of Multipath Rayleigh Fading: Multiuser 

In the multiuser case, we amukte the probalhlity of miss detection and probability of &lse 
alarm in a Rayleigh fading channel. There are three paths in the uplink. In this simulation we 
have assumed riming alignme nt [1] for all users which are under the power control of the base 
station. This essentially means that the first path of all users has i:ero or one chip-time relative 
delay. No inter-cell interference is considered. Power control is assumed and hence near-far 
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Probability of false aiarm & miss detection vs number of users in CFAR detection syste 



Figure 5.10; Probability of false alarm and miss detection versus the number of users 
in CFAR detection. 

effect due to cell-geometry is neglected. The assumptions are valid since the searcher carries 
out code-acqmsition for traffic channds. A MS can send traffic only -when other 
synchronisation as required in cell-search procedure has been achieved. Figure 5.10 shows the 
probability of miss detection for all the three paths and the total miss detection of the system. 
The simulation is carried out for 15 users. We observe that the probability of miss detection 
approaches a limit after a certain number of users become active in a system. This is due to the 
correlation of the Barker sequence part of the scrambling code. The output of the Barker CMF 
should be above the specified threshold and once this level is achieved, subsequent increase in 
signal power due to increased number of users having the same relative delay, does not change 
the probability of miss detection. In Figure 5.10 we observe that the false alarm rate is constant 
and is the same as that of the design parameter. 

5.5 Petfofmance of Code Acquisition by Tree Search 

In this section we analyse the performance of the tree search detection system which is used 
for user specific code acquisition. The signals arriving at the BS are Rayleigh faded and the 
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correlation of the scrambling codes due to larger number of users is considered Gaussian as 
per Central limit Theorem [45]. 

5.5.1 Inverse Filter 

The scrambling code is modified at the transmitter by decorrelating the 14* — 128* and 142"“^ — 
256* chip of the scrambling code with a decorrelator filter in order to reduce the false alarm in 
CFAR detector. The sequence will have to be restored back to its original form for tree search 
detection. This is achieved with the help of inverse filtering and defined in (4.23). In order to 
analyze the performance of the inverse filter, we sh^ compare the received waveform at the 
BS of the sequence defined by (4.8) and (4-15). The former sequence 




Figure 5.11: Comparison of inverse filter output signal with original sequence 

is without decorrelation and the latter sequence is with decorrelation and inverse filtered. The 
two sequences are obtained at the end of a multipath Rayleigh faded channel cormpted by 
AWGN. We also correlate the output of the inverse filter with the original sequence. We find 
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AWGN. We also correlate the output of the inverse filter with the n ngin al sequence. We find 
that the output of the inverse signal is highly correlated with the original sequence. The inverse 
filter restores back the decorrelated part of die scrambling sequence. 

5.5.2 Performance of MHT 

The output of the CFAR detector is used to enable a tree search block. The CFAR block also 
provides the value of which is the scaling factor for the threshold for MHT. In the present 
simulation, we have conadered multipath Rayld^ fading channel in the presence of Gaussian 
noise. The users have either two or three multipath with the first path either with 0 or 1 chip 
delay, second multipath of 2 or 3 dbip delays or third multipalh of 4 chip ddays. The design 
parameters chosen are = 0.05 , =0.1 and c, = Cj = 0.5 . We have used the parameters 

obtained from the partial corrdatioo bounds obtained in section 4.4.2. The plots for 
probability of miss detection for different mullipafhs for the number of active users in the 
system are shown in Figure 5.12. In Figure 5.13, we have the probability of false alarm for each 
of the multipath shown against variation of useis in the system. 

Comments: The value of the joint channel estimate obtained from the CFAR system is used 
to dynamically change the upper and lower threshold for the MHT. We have the following 
observations: 

• The relative mean power of the first path is -3 dB, the second path is -6 dB and the 
thkd path is -9 dB. We compare the equations for the threshokk in (3.61), (3.62) and 
(3.63) with that of the equation of a strai^t hne y = mK + c, where is the slope and c 
is the intercept We find that the intercept is divided by die factor y. For stronger paths 
the intercept will be small since y is higher for stronger paths. The two thresholds will 
be doser and steeper. This causes some of the test statistics to cross the upper 
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tiiteshold resuitkig in fake aknn* Tlius tkc fijtst path has a moderate miss detection 
probability. 

The third path has a converse approach than the first path. When there are fewer users 
in the system, the value of y is low and the thresholds are wider apart. Due to 
interference from stronger path signal, the test statistics cross the upper threshold. This 
results in the path having least probability of miss detection. With the addition of more 
users in the system, the value of y becomes higher resulting in steeper and narrower 
thresholds. This prevents the test statistics from crossing the upper threshold. Instead 
they clear the lower threshold and have higher probability of miss detection. This also 
results in lesser false alarm. 

• In comparison to the fiirst and the last delay paths, die second path shows a mixed 
performance. The dynamic thresholds hold good for most conditions and it 
approaches the design probabilities of false alarm and miss detection. 



Figure 5.12: Plot of Probability of miss detection versus number of users in tree 
search 
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Figure 5.13: Plot of Probability of false alarm versus number of users in tree search 
5.5.3 Average Sample Number 

The average sample number is the number of samples required by the TSPRT to make a 
decision. This factor depends on die values of r,and We have chosen die values to be such 
that the test is an equal mixture of FSS and SPRT. The plot is shown in Figure 5.14. 

Comments: The plot shows that when there are fewer users in the system, for padas with 

high er relative mean power, the threshold bang hi^er requires more number of samples while 
it is the other way round for weaker paths. As the interference rises, the weaka: path requires 
more samples to make a decision and the stcon^r paths require lesser number of samples- This 
is commensumbe with the results in Figure 5.14. Thus the number of samples tequiced using 
TSPRT is much lesser than that of a full search and the code acquisition is fast 
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Average sample number vs number of users in Tree Search 



Figure 5.14: Average number of samples required for detection 

Conclusioiis 

In this chapter we have seen that by concatenating the existing scrambling sequence with 
Barker sequence, we are able to detect the existing multipaths at tihe BS in 13 chips duration. 
By decorrelating the scrambling sequence between 14* — 128* and 142"^ — 256* chip with the 
inverse of Barker sequence we are able to contain die false alarm rate to the designed 
parameter. User specific code acquisition is achieved by tree search using MHT. The sequential 
test has reduced computational complexity and it requires lesser number of samples as 
compared to FSS. Therefore the acquisition is faster. If the test prolongs to the truncation 
stage as would be the case for weaker paths, then the total acquisition time required is 128 
chips. But for the stronger paths the acquisition is achieved in lesser number of samples. 
However, the pnce paid is that the probability of false alarm and nniss detection are higher. 
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Chapter 6 

Conclusions and Future Work 


This chapter summaiizes tihe tnajoi: contctbutioiis ia this thesis and presents possible fiitoce 
directions which coiiid be extension of the research work in this thesis. 

6.1 Conclusion 

It had been proved by Abraham Wald that sequential tests need mudh less sample size than a 
FSS test. Thus, systems employing sequential tests have a lower computational complexity than 
FSS tests. TSPRT requires much lesser sample size than a FSS for the same probabilities of 
false alarm and miss deteettoru We have proposed a new algorithm using the excellent 
autocorrelation properties of the Barker sequence along with tree search detection for the code 
acquisition problem. 

The probability of false alarm for multipath detection can be a chosen design parameter. For 
the tree search using MHT, the error probabilities obtained are less than or equal to that of the 
design parameters. The average sample size is for decision is much lower than the FSS test 
sample size. 

6.2 Future Work 

• For the CFAR detectioii system we have analysed the stahstical model for a Rayleigh 
Padmg rFannel. Si milar analysts needs to be done for Rician foding and log-normal 
shadowing channels. 
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• For MHT we have analysed imdef Gaussian channel conditions and for Rayleigh faded 
signals at the output of a CMF. It is suggested that sequential CFAE. detection be 
employed for tree search [55]. 

• The thresholds for TSPRT have been modified from the threshold derived for 
Gaussian distribution. It will be worthwhile to derive these for other distribution using 
asymptotic analysis. 

• The theoretical analysis for the acquisition time using TSPRT for paths with varying 
SNR can be made as gjiven in [56] [51]. 

• The existing system can be extended to multiple antenna systems using 2-D Rake 


receiver. 
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